On temporal variations of the multi-TeV cosmic ray anisotropy 
using the Tibet III Air Shower Array 
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ABSTRACT 

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV 
cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November 
to 2008 December. To explore temporal variations of the anisotropy, the data 
set is divided into nine intervals, each in a time span of about one year. The 
sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to 
year over the entire observation period of nine years. This indicates that the 
anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities 
since the observation period covers more than a half of the 23rd solar cycle. 

Subject headings: cosmic rays — diffusion — ISM: magnetic fields — solar 
neighborhood — Sun: activity 
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Introduction 



Galactic cosmic rays (GCRs) are high-energy nuclei (most protons) which are believed 
to be accelerated by supernova remnants (SNRs) in our Galaxy and continuously reach 
the Earth after propagating in the Galaxy and heliosphere. The intensity of GCRs 
is nearly isotropic due to deflections in the Galactic magnetic field (GMF). However, 
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The large-scale sidereal cosmic ray (CR) anisotropy may arise from several causes. 
Firstly, it may result from the uneven distribution of CR sources such as SNRs and the 
process of CR propagation in the Galaxy. Using the diffusion model, it is expected that 
the amplitude of anisotropy satisfies A ~ D(R) ~ i? 1 / 3 , where D{R ) is the diffusion 
coefficient which de pends on the rigidity R of charged particles (e.g. 



Shibata et al 



Berezinsky 



1990 



20041 ). Thus one can expect that the anisotropy amplitude will increase 
with energy. However, the observation of the anisotropy above hundreds of TeV is smaller 
than what is expected. It has been suggested that a nearby (within 1-2 kpc from the 



Sun) source in the opposite direction agai nst the globa 



expected large amplitude to some degree (jErlykin et al 



CR diffusion may suppress the 



20061 ). Other factors such as the 



complicated composition of CRs and is otropic CRs from the Galactic ha 



helpful for reducing the anisotropy (e.g. 



Erlykin et al. 



2006 



Ptuskin et al 



o can also be 



20061 ). Secondly, 



the anisotropy can also be induced through both large-scale and local magnetic field 
configurations, possibly including effects of the heliosphere. Therefore, it is a useful tool to 



probe the local i nterstellar space surrounding the heliosp 



the heliosphere (j Amenomori et al. 



2007 



lere and the magnetic structu re of 



Munakata et al. 



2008 



Amenomori et al. 



20091 ). In 
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addition, an expected anisotropy is caused by the relative mot ion between the observer a nd 



the CR plasma, known as the Compton-Getting (CG) effect (jCompton &: Getting 



By an alyzing the events recorded by the Tibet III Air Shower experiment, 



1935). 



Amenomori et al. 



( 120061 ) showed that GCRs corotate with the local GMF environment according to the null 



result of the Galactic CG effect. The CG modulation due to the Earth's orbital motion 



around the Sun has been successful 
energy ranges (e.g. 



Cutler & Groom 



detected by several GCR experiments in multi-TeV 



1986 



Amenomori et al 



2004 . 



2006 



Abdo et al. 
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From the analysis of numerous experiments, it can be seen that both the amplitude 
and the phase of the best- fit first harmonic vary with CR e nergy in a wide range from 



tens of GeV to PeV (IGuillian et al. 



2007 



Abbasi et al. 



20091 ). Below several tens of GeV, 



solar modulation effects are most notable for GCRs. GCRs interact with the solar wind 
magnetic field, both an ordered field and irregular field components, after entering the 
heliosphere. The spatial distribution of GCRs can reflect the magnetic structure in the 
solar wind. With increasing energy, CRs become less sensitive to the solar modulation. 
It is well known that the flux of GCRs with energy per nucleon in the energy range of 
~ 10 11 — 10 14 eV has a sidereal anisotropy of the order O(10 -3 ). The gyro radius of CRs 
in this energy range in a GMF of 3 /iG is about several AU— 0.03 pc, which is much smaller 
than the size of the Galaxy. In the multi-TeV range, the gyro radius of hundreds of AU 
becomes comparable to the spatial scale of the heli osphere in the nose direct ion toward the 



Washimi Tanaka 



J6) 



1996). However, it 



upstream side of the interstellar medium flow (e.g. 
is known that the heliosphere has a long heliotail, the modulation in the heliotail remains 
possible. Therefore, the large-scale sidereal anisotropy of CRs in this energy range gives 
us an important clue about the magnetic field structure of the heliosphere or the local 
interstellar space surrounding the heliosphere. 



The solar cycle shows a quasi-period of about 11 yr and the global magnetic polarity 
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reverses with a quasi-period of two solar activity cycles (e.g. ISolanki et al.l 120061 ). Since 
GCRs are modulated by solar activities in the given energy range mentioned above, it might 
be expected that the sidereal anisotropy may follow the variation of solar cycle. There 
have been many experiments devote d to the study of temporal variations of GCR sidereal 



anisotropy. 



Nagashima et al 



(1989) found that both the amplitude and the phase of the 
sidereal anisotropy with the rigidity ~ 10 13 V had no significant changes for more than ten 
years, except that the phase changed at the epoch (1979-1980) when th e polar magnetic 



field o f the Sun reversed its polarity. At primary energy as low as 2 TeV, 



Kozyarivsky et al. 



(119951 ) demonstrated that the sidereal anisotropy was constant within the accuracy of the 
experiment, with observations in the subsequent ten years (1982-1991) after this epoch 
of polarity reversal. Recently, underground muon observations showed that yearly mean 
harmonic vectors of the sidereal anisotropy in the sub- TeV region are more or less stable in 
phase and amplitude over 20 years since 198 5 without showing any significant correlation 



with the solar activity and magnetic-cycles (IMunakata et al. 



20081 ). However, Milagro 



experiment recently reported an amplitude increase of the sidereal anisotropy at 6 TeV in 
the lat ter half of the 23rd solar cycle (from 2000 July to 2007 July), while the phase remains 



stable ( lAbdo et al. 



20091 ) . The contradiction needs further checks from other experiments. 



Covering almost the same field of view and with similar sensitive energy range, the 
Tibet AS7 experiment seemingly observed no variation in the sidereal anisotropy using 
data from both Tibet H P and Tibet HI array w hich was divided into two parts by the solar 



maxim um around 2001 (lAmenomori et al. 



20061 ) . In the multi- TeV region, 



Amenomori et al. 



(l2005al ) found that both the amplitude and phase of the first harmonic vector of the daily 
variation are remarkably independent of primary energy. Taking advantage of the large 
field of view and high count rates as well as the good angular resolution of the incident 
direction, the Tibet III Air Shower Array provides currently the world's highest precision 
two-dimensional measurement of GCR intensity in this energy range. The observation 
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period runs from 1999 November to 2008 December, covering more than a half of the 23rd 
solar activity cycle from the maximum to the minimum. Therefore, we can do more precise 
study on the sidereal anisotropy variation year by year in correlation with the solar activity 
cycle. In this paper, we analyze temporal variations of sidereal anisotropy of multi-TeV 
GCR intensity using the data of Tibet III array from 1999 November to 2008 December. 
At this point, we cannot investigate the influence of the polarity reversal of the global solar 
magnetic field on the sidereal anisotropy due to the lack of data before the current magnetic 
field reversal. Note that we do not use the data of Tibet HD to perform the year-by-year 
analysis due to the limitation of the statistics, although it covers the period before the 
reversal. 



2. Tibet Air Shower Array Experiment 

The Tibet Air Shower Array experiment has been operating successfully at Yangbajing 
(90.522° E, 30.102° N; 4300 m above the sea level) in Tibet, China since 1990. The array has 
been gradually upgraded by increasing the number of counters from the Tibet-I array with 



65 p 



astic scintillation detectors placed on a lattice with a 15 m s pacing (Amenomori et al 



2002|). 



19921 ). The Tibet III array was completed in the late fall of 1999 (lAmenomori et al.l 
The array is composed of 497 fast timing (FT) detectors and 36 density (D) detectors, 
covering a surface area of 22,050 m 2 . Each FT detector equipped with a plastic scintillator 
plate and a 2 inch photomultiplier tube has a cross sectional area of 0.5 m 2 and is deployed 
at a lattice with a 7.5 m spacing. A 0.5 cm thick lead plate is placed on top of each counter 
in order to increase the detector array sensitivity by converting 7 rays into electron-positron 
pairs. A CR event trigger signal is issued when any fourfold coincidence occurs in the FT 
counters recording more than 0.6 particles, resulting in a trigger rate of about 680 Hz at 
a few-TeV threshold energy. The shower size pft is regarded as an estimator for the 
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primary particle energy, where the size of pft is defined as the sum of particles per 
m 2 for each FT detector. During 2002 and 2003, the inside area of the Tibet-Ill array 
was further enlarged to 36,900 m 2 by installing additional 256 detectors. This full Tibet 
III array has been operating successfully since 2003. In the present analysis, to keep the 
uniformity of the data quality throughout the entire observation period from 1999 to 2008, 
we reconstructed CR air shower events obtained using the configuration of detectors which 
was completed in the late fall 1999 even for the full Tibet-Ill array. 

In the present analysis, CR events are selected based on the following four criteria: (1) 
estimated air shower core location should be inside the array; (2) the zenith angle of the 
incident direction should be less than 45°; (3) any fourfold coincidence in the FT counters 
should record a signal of more than 0.8 particles; (4) when 10 < Yl Pft < 178 is satisfied, 
corresponding to a modal energy of about 5 TeV. 

In total, about 4.91 x 10 10 CR events are used in the present analysis. 

3. Analysis and Results 

Sitting on an almost horizontal plane, the Tibet III Air Shower Array has almost 
azimuth-independent efficiency in receiving GCR shower events for any given zenith angle. 
Therefore, the equi-zenith angle method is adopted in the present analysis. In brief, for a 
candidate "on-source window", GCR air shower events recorded simultaneously in the side 
band of the same zenith angle belt can be used to construct the "off-source windows" and to 
estimate the background events. This method can eliminate various detection effects caused 
by instrumental and environmental variations, such as changes in pressure and temperature 
which are hard to be controlled and tend to introduce systematic errors in measurement. 

With the large statistics, we can construct a two-dimensional sky map to reveal detailed 
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structural information of the large-scale GCR variations beyond the simple one-dimensional 
profile along the time scale. The idea of this method is that for each short time step (e.g. 8 
minutes in the present analysis), for all directions, if we scale down (or up) the number of 
observed events by dividing them according to their relative CR intensity, then statistically, 
those scaled observed numbers of events in a zenith angle belt should be equal anywhere. 
A total x 2 can be built accordingly, and the relative intensity of CRs I and its err o r Al in 



Amenomori et al 



fl2005bf ) 



each direction can be solved by minimizing the \ 2 function (see 
for more details). Given the fact that the detector can quickly scan the sky along right 
ascension direction as a result of the self rotation of the Earth, the intensities of CRs 
along right ascension direction are measured under the same condition for each declination 
belt. Therefore, the modulation along right ascension direction in each declination belt is 
accurate and irrelevant to the inhomogeneous efficiency of the detector. However, as having 



been pointed out by (e.g. 



Andreyev et al. 



2008 



Kozyarivsky et al 



20081 ). when comparing 



CR intensities among different declination belts, we need absolute efficiency calibration 
along declination direction which is unfortunately not possible for current experiments. 
Once we have the absolute CR intensity in declination direction Idec,iD{5), we will obtain 
the true two-dimensional anisotropy map, Itrue,2Di.oi,5) = I qU asi,2D(®,8) x Idec,iD(8), here a 
and 5 are right ascension and declination in the celestial coordinate respectively. Strictly 
speaking, the present analysis gives a combination of one-dimensional modulation curves in 
a fine set of declination bins - I qU asi,2D((x, 8) (For simplicity, we will omit the subscript of 
I qU asi,2D in the following text). 



1976 



Scargle 



19821 ) with 



Using Lomb-Scargle Fourier transformation method (jLombl 
CR data recorded by the Tibet III array, Tibet AS7 experiment showed that besides the 
well-known solar diurnal, sidereal diurnal and sidereal semi-diurnal modulations at a level 
of ~ 10 -3 , no other periodicity was found to have high en ough significan ce from 1 hour to 2 



years in the energy range from ~ 3.0 TeV to ~ 12.0 TeV (ILi et al. 



20081 ). The sidereal daily 
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variation can be described by the first and second harmonics and the solar daily variation 
can be described by the first harmonic alone. When applying the above mentioned method 
to data, the fit function is natural to contain both sidereal time and solar time modulation 
components. The other alternative method we used to separate modulations was to properly 
fold the data according to the sidereal time or solar time periodicity. This method was 



used in (lAmenomori et al. 



2005b 



20061 ) when we analyzed CR intensity variations in the 
sidereal time frame with several years data samples. Using this method, each event had to 
be re-weighted to form an exactly one year long and uniformly distributed data series before 
folding. Therefore, the disadvantage of this method is that more than one year's data are 
needed, which can not be followed in current work. In this case, we adopt a new method 
of fitting all modulation components simultaneously to perform the year-by-year analysis of 
CR anisotropy. 

Based on the results of Lomb-Scargle Fourier transformation, in the present analysis, 
we assume that at any moment t, the relative intensity of CRs at any given direction (9, (ft) 
in the horizontal coordinate, is modulated as a product of I S id{<^sid,S S id) and I so i(oi so i, S so i). 
Here (a sic i, 5 sic i) and (a so z, <5 so z) are positions corresponding to the celestial coordinate in 
the local sidereal time frame and the local solar time frame of the same point (t, 6, 0) in 
horizontal coordinate, I S id and I so i denote the CR intensity in the sidereal time frame and 
the solar time frame respectively. So substituting / by I S id( a sid,$sid) x Isoii^soi^soi), the 
total x 2 can be written as 



x 2 = E 



2 

NobsMS. E<^<// [^obs(t,e^')/(i sld (< ld ^ sld )><isoiKoi' s 'soi))} 

t,8,4> \ "I I sid (a sid ,5 sid )xI 30 i(a ao i,8 sol ) Z)<^0' 1 

r 



v J N obs (t,8,& | S^^K 63 (f,e,0 , )/(/ s 2 Id K Id ,.5L: d )x J L« o P- ; s Oi ))] i j m 

\^ d (^ id ,S aid )^ ol (a aol ,S a oi) + (E^'l) 2 f I ' 11 

Here, N obs (t,9,(fi) denotes the number of observed events in the "window" (9,4>) in the 
horizontal coordinate at the moment t. After minimizing the \ 2 function, I S id(a S id, 5 S id) and 
hoi^soh^soi) can be obtained simultaneously. 
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Fig. 1. — Sidereal diurnal variation of CR relative intensity with the representative primary 
energy of 5 TeV averaged over all nine phases of Tibet III Air Shower Array from Nov. 1999 



to Dec. 2008. (A) 2D maps obtained with the present method mentioned as 



(B) 2D maps obtained with the former method used in (lAmenomori et al.ll2006l ); (C) the ID 



Equation ([!]); 



projection of 2D maps averaged over all declinations for comparison. 

The data were acquired by the Tibet III Air Shower Array for 1915.5 live days from 
Nov. 1999 to Dec. 2008. It covers all nine running phases of the Tibet III Air Shower 
Array. As described in Equation (TjQ), for any direction in the horizontal coordinate at any 
moment in the observation period, a x 2 function can be constructed. Covering all nine 
phases of Tibet III array, we got a total \ 2 accumulated over the entire observation period. 
By minimizing it, the CR intensity maps in the frame of the local sidereal time and the 
local solar time averaged over all nine phases were obtained simultaneously. Fig. [TJ (A) 
shows the intensity map of GCRs with the modal energy of 5 TeV in the local sidereal 



13 



time frame averaged over the whole observation period of Tibet III array. The result is 
consistent with former observation re sults of Tibet III experiment by different methods 



( [Amenomori et al. 



2005a . 



2006 



20071 ). The amplitude of the sidereal anisotropy is about 



0.1%, and the maximum around 6 hr in the local sidereal time. 



"loss-cone" anisotropy components (e.g. 



Nagashima fc Mori 



1976 



he so-called "tail-in" and 



Nagashima et al 



1998|) 



are clearly shown in the map. The maximum phase shifts to earlier hours as the viewing 
declination moves southward in two-dimensional intensity map. The excess around the 
Cygnus arm direction can also be seen in the map. 

For a detailed examin ation between the present method and the former method used 



in 



Amenomori et al. 



( 120061 ). we constructed a two-dim ensional map of CR relative intensity 



with former method used in 



Amenomori et al. 



( 120061 ) with the same data samples, showed 



in Fig. [JJ (B). The ID projections of Fig. [TJ (A) and (B) averaged over all declination belts 
(—14.89° to 75.11°) are displayed in Fig. [1] (C). From the comparison, it can be seen that 
the results agree with each other well. 

To study the temporal variation of the sidereal anisotropy, the data was divided into 
nine subsets, corresponding to nine running phases of Tibet III array each in a time scale 
of about one year, as summarized in Table [TJ 

Using data samples recorded during each separated phase, the GCR relative intensity 
maps in the local sidereal time frame averaged over each phase were obtained according 
to Equation (pQ). The two-dimensional (2D) relative intensity maps of GCRs with modal 
energy around 5 TeV in the local sidereal time frame of different phases are displayed in 
Fig. [2j and the corresponding one- dimensional (ID) projections over all declinations are 
also shown. The solid red markers in each plot denote the relative intensities of GCRs 
in the local sidereal time frame over the corresponding observation period, while the blue 
dashed smooth curves in plots represent variations of GCR relative intensity averaged over 
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Table 1: Definition of nine phases of Tibet III from 1999 November to 2008 December 
Phase Start time End time Live days Number of used CR events 



1 


Nov. 


18, 


1999 


Jun. 


29, 


2000 


173.1 


5.16x10 


2 


Oct. 


28, 


2000 


Oct. 


11, 2001 


283.7 


8.14x10 


3 


Dec. 


05, 


2001 


Sep. 


19, 


2002 


201.8 


5.59x10 


4 


Nov. 


18, 


2002 


Nov. 


18, 


2003 


259.1 


6.34x10 


5 


Dec. 


14, 2003 


Oct. 


10, 


2004 


123.6 


3.07x10 


6 


Oct. 


19, 


2004 


Nov. 


15, 


2005 


277.6 


6.79x10 


7 


Dec. 


07, 


2005 


Nov. 


03, 


2006 


114.5 


2.71x10 


8 


Nov. 


06, 


2006 


Feb. 


25, 


2008 


269.2 


6.36x10 


9 


Mar. 


02, 


2008 


Dec. 


03, 


2008 


212.9 


4.91x10 



all nine phases of Tibet III array, the same as the bottom plot with the marker of "present 
method" in Fig. [TJ From the comparison of GCR sidereal anisotropy in different phases 
from 1999 November to 2008 December, it can be seen that the CR intensity variation in 
the local sidereal time appears fairly stable year by year. 

Furthermore, we take a \ 2 test to check the consistency among different phases. As 
shown in Equation ([TJ), the CR intensity variation along the local sidereal time averaged 
over all nine phases contains contributions of each single phase. To avoid the correlation 
between the two plots which are used in the test, we compare the result of each single 
phase with the average one of all the other eight phases except the single one instead 
of the average one of all nine phases. The obtained x 2 / n df value and the probability of 
each comparison were labeled in Fig. |2j The test also indicates stability for the sidereal 
anisotropy with time. 

The observation period of Tibet III Array covers more than a half of the 23rd solar 
activity cycle from the maximum to the minimum. So it implies that the sidereal anisotropy 
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Fig. 2. — Cosmic ray intensity variation in the local sidereal time frame for CRs with the 
modal energy around 5 TeV in the 9 phases of Tibet III Array. Top: 2D intensity map of 
each phase; Bottom: ID projection averaged over all declinations. In bottom plots of each 
panel, the red crosses in each plot show the intensity variation over each phase respectively, 
while the dashed blue lines represent the intensity averaged over all nine phases of Tibet III 
array 



of multi-TeV GCRs i s insensitive to th e solar activity It disagrees with the recent result of 



Milagro experiment (lAbdo et al.ll2009l ). which shows an increase in the amplitude of the 



sidereal anisotropy with time while the phase remains stable. 
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4. Conclusions 

In this work, we investigate temporal variations of the large-scale sidereal anisotropy 
of GCR intensity using the data of Tibet III Air Shower Array from 1999 November to 
2008 December. Totally ~ 4.91 x 10 10 CR events are used. The data is divided into nine 
intervals, each in a time span of about one year. We find that, in the multi-TeV energy 
range, the sidereal anisotropy is fairly stable year by year over all nine phases of Tibet III 
Array, which covers more than a half of the 23rd solar cycle from the maximum to the 
minimum. It indicates that the anisotropy in this energy range appears insensitive to solar 
activities. This feature can give some constraints on the origin of the sidereal anisotropy, 
which has no convincing and widely accepted explanations so far. 
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